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Abstract 

The analytic properties of Nilsson's Modified Oscillator (MO), which was first introduced 
in nuclear structure, and of the recently introduced, based on quantum algebraic techniques, 
3-dimensional g-deformed harmonic oscillator (3-dim g-HO) with Ug(3) D sOq(3) symmetry, 
which is known to reproduce correctly in terms of only one parameter the magic numbers of 
alkali clusters up to 1500 (the expected limit of validity for theories based on the filling of 
electronic shells), are considered. Exact expressions for the total energy of closed shells are 
determined and compared among them. Furthermore, the systematics of the appearance 
of supershells in the spectra of the two oscillators is considered, showing that the 3-dim 
g-HO correctly predicts the first supershell closure in alkali clusters without use of any extra 
parameter. 
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1. Introduction 

Supershells, which are seen as beats of the deviation of the total energy of many particle 
systems from the function describing its average behavior vs. the number of particles N, 
are known to be a general property of the spectrum of potentials having sharp edges [^. 
Supershells in metal clusters have first been studied by Nishioka, Hansen and Mottelson 
in terms of phenomenological mean field potentials. 

In the other hand, using recently developed quantum algebraic techniques 0, it has 
been shown |^ that the magic numbers appearing in alkali clusters can be successfully 
reproduced up to 1500 (which is the expected limit of validity of theories based on the 
filling of electronic shells 0) by the 3-dimensional g-deformed harmonic oscillator (3-dim 
g-HO), which possesses the Ug(3) D sOg(3) symmetry Furthermore, the magic numbers 
appearing in several divalent (Zn, Cd) and trivalent (Al, In) metal clusters have been 
satisfactorily reproduced @] by the same model, in terms of only one free parameter, the 
deformation parameter r (with q = e^, where r is a real number). It is therefore of interest 
to examine if the 3-dim g-HO can predict supershells and which these predictions are. It 
should be noticed that the calculation of supershells in the framework of the 3-dim g-HO 
will be parameter free, since the single parameter of the model has been fixed in reproducing 
the magic numbers for each kind of clusters 0. 

In addition to the 3-dim g-HO, Nilsson's Modified Oscillator (MO) [0, §, which has 
first been used in describing the structure of atomic nuclei, has also been early employed 
in describing atomic clusters ||^ (after dropping the spin-orbit interaction, which plays an 
essential role in nuclear structure but is absent in the case of atomic clusters). It is therefore 
of interest to study the possible appearance of supershells in the framework of this model 
as well. 

For the determination of supershells the method of Ref. can be employed. Before 
doing so, one has however to examine the analytic properties of the spectra of the two 
oscillators, in order to be able to apply meaningful truncation schemes. Furthermore, the 
average behavior of the total energy of a system of many particles (an atomic cluster in 
the present case) as a function of the particle number N should be determined, since it is 
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needed in the procedure of the study of supershells. As a result, these tasks will be carried 
out for both models, before any attempt for the determination of supershells is made. 

Supershells have been early predicted by Balian and Bloch |jTO| in the study of electrons 
moving in a spherical cavity, which by analogy can be used for the valence electrons in 
a metal cluster. The comparison of the stringent predictions of the theory of Balian and 
Bloch for various characteristics of the supershells to the results of the present models 
turns out to be a fruitful testing procedure. 

In Section 2 the analytic properties of Nilsson's Modified Oscillator will be considered, 
while the corresponding properties of the 3-dim g-HO will be studied in Section 3. In 
Section 4 supershells in Nilsson's Modified Oscillator will be studied, while supershells in 
the framework of the 3-dim g-HO will be considered in Section 5. Finally in Section 6 a 
discussion of the present results and plans for further work will be given. 

2. Nilsson's Modified Oscillator (MO) 

The potential of the Modified Oscillator (MO) introduced in nuclear physics by Nilsson 
0, ^ (with the spin-orbit term omitted) is 



V=hu;p'-hu;fi'{L^-<L^>r,), p = r^^, (1) 

where 

< V >„= (2) 

The effect of the term is to flatten the bottom of the potential. In addition it causes an 
overall compression of the spectrum, which is avoided through the addition of the < 
term, which will be further discussed later. 

The energy eigenvalues of Nilsson's modified harmonic oscillator are 0, ^] 

Enl = nuj(^n + ^^- hufi' + 1) - ^n{n + 3)) , (3) 

where n is the number of vibrational quanta and / is the eigenvalue of the angular momen- 
tum, obtaining the values / = n, n — 2, . . . , or 1 (depending on n being even or odd, 
respectively) . 
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The number of states having energy (n + |^ hu! in the case of even n is 

± p,.,).(!l±lK!l±^. (4) 

l=0,l=even 

where only the even values of I are included in the summation and the sum 

t^^"-^ (5) 

X=l ^ 

has been used. The same result is obtained for odd n, in which case only the odd values of 
I are included in the summation. Taking into account the spin of the particles the number 
of states having energy [n + |^ hu is 

712= (n + l)(n + 2). (6) 

The sum of the eigenvalues of \? within each shell in the case of even n is 

L,= t i(; + l)(2i + l) = "<"+^<" + ^"" + ^'. (7) 

i=0,i=even 

where only the even values of / have been included in the summation and, in addition to 
Eq. (5), the sums 

^^,_ n(n + l)(2n + l) 

and 

(9) 

have been taken into account. The same result is obtained for n being odd. Taking the 
spin of the particles into account one has 

L, = 2L. = "'"+'>'" + '"" + ^^ (10) 

Thus the average per particle of the square of the angular momentum within each shell is 

< >„= £l = £a = !!(!L±1), (11) 

Hi 77-2 2 

a result which has already been used in Eqs. (1) and (3), in order to keep the "center of 
mass" of each shell constant, i.e. to counterbalance the overall compression of the spectrum 
caused by the term alone. 



The total number of particles which can be accommodated in the levels of the shells up 
to the n-th shell included, is 

^^i:(^+i)(.+2)= '"+^)'"+^)("^^) . (12) 

x=0 

where the spin of the particles has been taken into account and Eqs. (5), (6), (8) have been 
used. It should be remembered throughout the present work that N stands for the total 
number of particles, while n stands for the number of vibrational quanta. 

The contribution of the first term of Eq. (3) to the total energy of the particles up to 
the n-th shell included (and taking the spin of the particles into account) is 

E(n) = B., t + ^) + IX'^ + ^) = t ^ " • (") 

where Eq. (6) is used for the degeneracy within each shell and Eqs. (5), (8), (9) have been 
used for performing the summations. For later use we notice that omitting the ground state 
energy contribution in a similar manner one finds 

X ^ A / X/ ^^ * n(n+l)(n + 2)(n + 3) 
E'{n) ^hu;Y,x{x + l){x + 2) ^hu^ — ^, (14) 

x=Q ^ 

while a perturbation in the energy would have given an additional term 

/ X ^- 2/ .X/ nx n(n+ l)(n + 2)(n + 3)(4n+ 1) ,^ , 

E2{n) ^?iuj^x^{x + l){x + 2) ^huj^ ^ ^, (15) 

x=0 

where the sum 

" 4 n(n + l)(2n + l)(3n2 + 3n-l) 

= ^ (16) 

has been used in addition to Eqs. (5), (6), (8), (9). 

The contribution of the second term of Eq. (3) to the total energy of the particles up 
to the n-th shell included is found by using Eq. (10) 

. " a;(a; + l)(a; + 2)(a; + 3) 

x=0 ^ 

while the contribution of the third term of Eq. (3) will be 

1 " 

E^{n) = +-nujii' E x{x + 3)(a; + l)(a; + 2), (18) 

2 x=0 



where use of Eq. (6) has been made. We remark that E^{n) and E/^{n) cancel. Thus we 
conclude that in Nilsson's MO only the first term of Eq. (3) contributes to the total energy 
of the particles up to the n-th shell included. 

The average energy per particle (up the n-th shell included) is then found using Eqs. 
(12) and (13) to be 

Ein) ^ /3 3\ 
< E >= —-^ = huj[-n + -] , (19) 



V4 2/' 

i.e. the average energy per particle is increasing linearly with the shell number, which is 
the number of vibrational quanta, as it is expected for a harmonic oscillator, since the 
angular momentum terms make no contribution, as we have already seen. The same result 
is obtained from Eq. (14) 

E'in) ^3 

< E' >= — ^ = hu-n, 20 

where the 3/2 term has been omitted already in Eq. (14), while a perturbation in energy, 
as seen from Eqs. (12) and (15), would have given an extra term 

< E^in) >= = hcu^n{An + 1), (21) 

which naturally contains a term. 

The lower part of the spectrum of Nilsson's Modified Oscillator, calculated from Eq. (3), 
is shown in Fig. 1( function of the parameter fi', together with the magic numbers 

appearing at different parameter values. The following observations can be made: 

a) The magic numbers at the left end of the figure are the ones of the 3-dimensional 
isotropic harmonic oscillator. 

b) The magic numbers appearing around /i' = 0.04 are in agreement with the magic 
numbers appearing in alkali clusters, up to 138 (see Ref. @] for more details). The agree- 
ment is destroyed beyond this point. 

c) Around the parameter value fi' = 0.02 the magic numbers up to 138 are a mixture 
of magic numbers of the 3-dimensional isotropic harmonic oscillator and magic numbers 
appearing around ^ = 0.04 (magic numbers of alkali clusters). 

d) Eq. (3) can be rewritten as 

Eni = ^hio + huon il + ^ j + hujn^^ - huj^'l{l + 1), (22) 



which clearly shows that in the case of /i' > the levels within each oscillator shell (char- 
acterized by a given value of n) are ordered according to the value of /, with the levels with 
higher values of / lying lower in energy, because of the last term in Eq. (22). This is indeed 
the case in Fig. 1(a), although the levels have not been labelled by the quantum numbers 
n, / because of lack of space. 

e) The level with / = n in particular lies lowest in energy within each shell and in general 
its energy is decreasing with increasing /i', since in this case Eq. (3) takes the form 

Enn = ^^00 + hujn ^1 + ~ Tiuv?^, (23) 
its derivative with respect to /i' being 

which is indeed negative for n > 1. Indeed the levels which lie lowest within each oscillator 
shell in Fig. 1(a) are the levels with / = n, which also show negative slope with increasing 

f) The fact that the ri^ term in Eq. (23) appears with a negative sign (for /i' > 0) can 
cause difficulties if one tries to describe a system with a large number of particles in terms 
of this oscillator. The derivative of Enn with respect to n 

f!^.S.(l4-,'„) (25) 

remains positive for 

" < ^ + i. (26) 
Beyond this value of n the derivative is negative, meaning that levels with higher values 
of n will lie lower in energy, making it difficult to define a cut-off for the number of shells 
taken into account. For example, if /i' = 0.04 (a value which has been found relevant 
to the description of metal clusters), the derivative remains positive if < 25.5. It is then 
clear that a reasonable truncation of the spectrum is possible only if the number of shells 
to be taken into account is less than 26 (taking the n = shell into account), otherwise 
no truncation is possible. This drawback of Nilsson's MO does not have any consequences 



in the case of nuclear structure, where the model has been first introduced, because of the 
small number of particles involved there (for which including the shells up to n = 8, shown 
in Fig. 1(a), suffices), but it can cause difficulties if one tries to employ this model for the 
determination of supershells in metal clusters, as we shall see later in Sec. 4. 

g) As an extention of f), one sees from Eq. (23) that Enn remains positive if n < ^ + 1- 
Thus, in the case of n' = 0.04 one should have n < 51, otherwise energies lower than the 
ground state energy will occur. 

3. The 3-dimensional g-deformed harmonic oscillator (3-dim q-HO) 
The space of the 3-dimensional g-deformed harmonic oscillator consists of the com- 
pletely symmetric irreducible representations of the quantum algebra Ug(3). In this space 
a deformed angular momentum algebra, sOq(3), can be defined 0. The Hamiltonian of the 
3-dimensional g-deformed harmonic oscillator is defined so that it satisfies the following 
requirements: 

a) It is an sOg(3) scalar, i.e. the energy is simultaneously measurable with the g-deformed 
angular momentum related to the algebra sOq(3) and its z-projection. 

b) It conserves the number of bosons, in terms of which the quantum algebras Ug(3) and 
sOg(3) are realized. 

c) In the limit g — > 1 it is in agreement with the Hamiltonian of the usual 3-dimensional 
harmonic oscillator. 

It has been proved that the Hamiltonian of the 3-dimensional g-deformed harmonic 
oscillator satisfying the above requirements takes the form 

where J\f is the number operator and C^^-* is the second order Casimir operator of the 
algebra sOq(3), while 

M, = ^ (28) 
is the definition of g-numbers and g-operators. 

The energy eigenvalues of the 3-dimensional g-deformed harmonic oscillator are then 

E,M = /.u;o |[n],g"+^ - ^-^^^^lU + 1].} , (29) 



where n is the number of vibrational quanta and / is the eigenvalue of the angular momen- 
tum, obtaining the values / = ra, — 2, . . . , or 1. 

In the limit of g ^ 1 one obtains \img^iEg{n,l) = hujQn, which coincides with the 
classical result. 

For small values of the deformation parameter r (where q = e"^) one can expand Eq. 
(29) in powers of r obtaining 

Eq{n, I) = hujQTi — UujqT + 1) — n{n + 1)) 

- huoT^ + 1) - ^n{n + l){2n + 1)) + 0(r^). (30) 

The number of states characterized by a given value of n, i.e. the number of states in 
the n-th shell, is still given by Eq. (4) if the spin is not taken into account, and by Eq. (6) 
with spin taken into account. 

The total number of particles which can be accommodated in the levels of all shells up 
to the n-th shell included is still given by Eq. (12), with spin taken into account. 

The analogue of the sum of the eigenvalues of within each shell in the case of even / 
is, in analogy with Eq. (7), given by 

n -■ 

^i»= E [/y/ + 1U2/ + 1) = 

l=0,l=even W V J W H ) 

((2n + l)(g2"+^ + g-2"-^) - (2n + 5)(g2"+i + g-^^^^) + 5(g + q-') - (q' + g'^) 

- ^^^^{q' + q~' + q' + g"^ ' 2q - 2g-^)), (31) 

where [l]q[l + l]q are the eigenvalues of the second order Casimir operator of sOg(3), and 
use of (2/ + 1) has been made for the degeneracy within a shell without taking spin into 
account. In performing the relevant summations one needs, in addition to Eq. (5), the 
sums 

n pr(n+l) _ i 

Ee-= , (32) 

x=0 ^ ^ 

n 1 

y xe™ = -r f ne"("+2) _ ^ i)e-("+i) + e^) , (33) 

x=o (e" - 1)' ^ ^ 
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of which the first is a simple geometric series, while the second can be derived from the first 
through differentiation with respect to the parameter r. One can easily see that for odd I 
a result identical to the one given in Eq. (31) is obtained. 

Using the definition of the g-numbers given in Eq. (28) the above result can be rewritten 
in the form 



(g2 - 1) 



(2n + 1) 



n + 



-4 



n + 1 



n 



+ 1 (n + 1) 



+ 



(34) 



or equivalently 
LiM = — 



n 



- + l)(n + l)[2], + — 



where, for example, use of the identity 

[n\g2 



q2 _ q 2 



— 1 9 —9 

q — q q^ — q 



(35) 



(36) 



has been repeatedly made. 

Eq. (31) can be rewritten in yet another form by using the definition for Q-numbers 
(see, for example, the review article in for relevant details) 

g"- 1 



Q-l 



where 



Using this definition Eq. (31) takes the form 



(37) 



(3^ 



(g-i)(g2_i)2 



{[n]QQ'/' {{Q' - 5) + 2n{Q' - r 



2(g-i)(g2-i) 



-n]Qg-i/2((g-2-5) + 2n(g-2-i 

n^Q + lf-niQ-l)'' 



(39) 
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In the limit q ^ 1, keeping terms of order up to one can see that Eq. (34) is reduced 
to Eq. (7), i.e. it is in agreement with the non-deformed case. For this calculation one 
finds helpful the Taylor expansion of g-numbers P] 

[n]g = n±—{n-n^) + —{7n-10n'^ + 3n^)± — — {31n - 49n^ + 21n^ - Sn") + ..., (40) 

oDU XOX^U 

where the upper signs correspond to q being a phase factor {q = e*"^ with r being real), 
while the lower signs correspond to q being real (g = e"^ with r being real), as in the present 
case. 

One can now proceed to the calculation of the total energy of the particles up to the 
n-th shell included. Using the identity 

= QMq, (41) 

where g-numbers of Eq. (28) (Q-numbers of Eq. (37) ) are used in the left (right) hand 
side and Q = q^ (Eq. (38) ), one finds that the contribution of the first term of Eq. (29) to 
the total energy is 

n n 

E,^,{n) = /icuo 5][x],g"+^(x + l)(x + 2) = huo ^ Q[x]Qix + l)(x + 2) 

x=0 x=0 

= ^^oj^^rjy ((^ + + 2)Q"+' - 2{n + l){n + 3)Q"+2 + (n + 2)(n + 3)Q"+i - 2 

' ^^° 3(Q^-1) ^^ + 1)(^ + 2){n + 3), (42) 
where, in addition to Eqs. (5), (6), (8), (32), (33) one also needs to use the sum 

n 1 

.=0 (e--i)3V y 

(43) 

which is derived from Eq. (33) by differentiation with respect to the parameter r. Using 
Eq. (37) one can easily see that Eq. (42) can be put in the more symmetric form 

((n + l)(n + 2)[n + 3]q - 2{n + l)[n + 2]Q{n + 3) + [n + l]Q{n + 2){n + 3)) 
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- ^^° 3(g^-l) ^'' + l)(^ + 2)(n + 3). (44) 
In the limit of Q ^ 1, keeping terms up to (where Q = = = e^"^ and thus T = 2r) 
one finds that Eq. (44) agrees with Eq. (13) of the non-deformed case. In this calculation 
it is helpful to use the Taylor expansion of Q-numbers 

No = ^ + ^(^^ + - + 1) + —{n^ -2n^ + n^) + ... (45) 

The contribution of the second term of Eq. (29) to the total energy is found in a similar 
manner. One has 

EiAn) = -nuJo2 — 2^^Li,,(x) = -^^o ^g2 _ 1)3 

{Q^[n]Q{2{Q^ - l)n + {Q^ - 2Q - 7)) - Q^[-n]Q{2{Q-^ - 1)N + {Q-^ - 2Q-' - 7)) 
_ l.n\n + 6)(Q + 1){Q^ - 1)^ + ln(Q + 1){Q^ + QQ' + 3AQ^ + QQ + 1)), (46) 

D 

where Eqs. (5), (8), (31), (32), (33) have been used and the spin of the particles has been 
taken into account. 

The lower part of the spectrum of the 3-dimensional g-deformed harmonic oscillator, 
calculated from Eq. (29), is shown in Fig. 1(b) as a function of the parameter r, together 
with the magic numbers appearing at different parameter values, while in Fig. 1(c) the full 
spectrum up to about 1500 particles is exhibited. The following comments and comparisons 
to Nilsson's MO are now in place: 

a) The magic numbers at the left end of both figures are the ones of the 3-dimensional 
isotropic harmonic oscillator, as in Fig. 1(a). 

b) The magic numbers appearing around r = 0.04 are in agreement with the magic 
numbers appearing in alkali clusters, up to 1500 (see Ref. for more details), which is 
the expected limit of validity for theories based on the filling of electronic shells 0, while 
in the case of Nilsson's MO the agreement is limited to the magic numbers up to 138 (Fig. 
1(a) ). 

c) Around the parameter value r = 0.02 the magic numbers up to 138 in Fig. 1(b) are 
a mixture of magic numbers of the 3-dimensional isotropic harmonic oscillator and magic 
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numbers appearing around r = 0.04 (magic numbers of alkali clusters). Beyond 138 other 
magic numbers appear. 

d) Numerical calculations show that in the case of r > the levels within each oscillator 
shell (characterized by a given value of n) are ordered according to the value of I, with the 
levels with higher values of / lying lower in energy, because of the last term in Eq. (29). 
This is indeed the case in Figs. 1(b) and 1(c), although the levels have not been labelled 
by the quantum numbers n, I because of lack of space. The ordering of the levels within 
each shell is the same as the one appearing in the case of Nilsson's MO (Fig. 1(a) ). 

e) As in Nilsson's MO, the level with I — n lies lowest in energy within each shell. 
However, in the present case the energy of this level is not decreasing with increasing r. 
This is true for all levels of the 3-dim g-HO: Their energies increase with increasing r 
(except in the cases of the levels with (n, I) — (0, 0) and (1,1), the energies of which remain 
constant with increasing r). 

f) As a consequence of e), no difficulties related to truncation appear in the present 
case. Stopping the level scheme at the / = n level of a given shell and taking into account 
all levels with lower n (i.e. all levels of the shells lying below the given one), one makes sure 
that all levels up to the given level have been included. Therefore in the 3-dim q-RO rehable 
truncations can be made, allowing for the description of systems with many particles. This 
point is one of the main advantages of the 3-dim g-HO in comparison to Nilsson's MO. 

A few more comments on the comparison between the 3-dim q-HO and Nilsson's MO 
are also in place: 

a) The < >„ term in Nilsson's Hamiltonian (Eq. (1)) has been put in "by hand" 
in order to guarantee that the "center of mass" of each shell will remain constant, so that 
shells will not be compressed because of the presence of the term. In the case of the 
3-dim g-HO it is clear that the opposite effect is present: The shells are expanded, because 
of the extra terms added by the g-deformation. One way to see this is by comparing the 
first order corrections appearing in Eq. (30) to the last two terms in Eq. (3). In both cases 
the l{l + 1) term causes compression of the shells, while expansion of the shells is caused 
by the n{n + 1) term in the case of the 3-dim q-HO and by the n{n + 3)/2 term in the case 
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of Nilsson's MO. Since the difference of these terms is 

n{n + 1) — -^n{n + 3) = -n{n — 1), (47) 

which is positive for n > 1, it is clear that the expansion in the case of the 3-dim g-HO 
will be stronger than the expansion in the case of Nilsson's MO, which is exactly balanced 
by the compression caused by the /(/ + 1) term. As a result, in the case of the 3-dim 
g-HO net expansion of the shells will occur, which will be of first order in the parameter r. 
Comparison of Figs. 1(a) and 1(b) clearly shows this effect. 

b) In Nilsson's MO the last two terms in the energy expression of Eq. (3) give no net 
contribution to the total energy up to the n-th shell included, and therefore to the average 
energy per particle as well, which turns out to be proportional to n (see Eq. (19) ). In the 
3-dim g-HO this dependence of the average energy per particle on n is given by the lowest 
order contribution from Eq. (44), while the next order contribution from Eq. (44), as well 
as the lowest order contribution from Eq. (46), give terms with higher powers of n. This 
property will have to be taken into account when supershells will be considered. 

c) In both models we have derived exact expressions for the total energy up to the 
n-th shell included, i.e. for systems of particles filling complete shells. In order to consider 
systems of particles for which the last shell is not full, one has to consider numerical 



methods, as the Strutinsky method |ri|, [1^, which are beyond the scope of the present 
study. 

4. Supershells in Nilsson's Modified Oscillator 

For studying the existence and properties of supershells in Nilsson's MO we are going 
to use the procedure employed by Nishioka et al. 0. For a given number of particles N 
the single particle energies Ej{n, I) of the N occupied states are summed up 

E{N) = j2E,{n^). (48) 

3=1 

This sum is then divided into two parts: A smooth average part Eav and a shell part E shell, 
which will exhibit the supeshell structure 

E{N) = Ea.{N)+Esheii{N). (49) 
14 



For the average part of the total energy a Liquid Drop Model expansion is used 



(50) 



This expansion should be adequate in the case of Nilsson's MO, for which the average 
energy per particle increases linearly with (see Eq. (19) ), but it should not suffice in 
the case of the 3-dim g-HO, for which the average energy per particle contains higher order 
terms, as we have seen in the previous section. For the latter case an expansion going up 
to a A^^ term 



should be more appropriate. In order to keep the calculations uniform and thus facilitate the 
comparisons between the two oscillators, we opted for using the expansion of Eq. (51) in all 
cases, although for Nilsson's MO the first three terms would have been adequate. Therefore 
in the case of Nilsson's MO very small values will be expected for ae, the coefficient of the 
N'^ term. 

The results for Egheii obtained with Nilsson's MO for 8 different values of the parameter 
/i' (assuming huj = 1) are shown in Fig. 2, while the relevant parameter values and rms 
deviations are exhibited in Table 1. As far as the way of calculation is concerned, the 
following comments apply: 

a) In all cases the maximum shell used was rimax = 30, implying that the last level 
included in the calculation was the one with {n,l) = (30,30), in order to ensure that the 
complete spectrum up to this point has been taken into account, as discussed in Sees. 2 
and 3. 

b) From Eq. (26) one sees that for fi' = 0.02 one should keep n < 50.5, while for lower 
values of fi' the limiting value of n lies even higher. This means that including the shells 
up to n = 30 is a reasonable truncation. 

c) The procedure of the calculation was as follows: First the summations described 
by Eq. (48), resulting in the total energy E{N) for each particle number A^, have been 
performed. Subsequently, in order to reduce the size of the calculation approximately by 



(51) 
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a factor of 10, the average E{N) was calculated every 11 points (i.e. for N = 6, 17, 28, 
. . . ), up to the point immediately below the cut-off of (n, /) = (30, 30), which is reported in 
Table 1 as N^ax- These averaged values of E{N) were subsequently fitted by the expansion 
of Eq. (51), resulting in the determination of Eav{N) at these points. Finally Esheii{N) has 
been obtained at these points as the difference E(N) — Eav{N) and plotted in Fig. 2. 
On the contents of Table 1 the following comments can be made: 

a) The behavior of the parameters as functions of fi' is rather smooth, with the exception 
of oe, the coefficient of N"^, which assumes very small values, as expected from the comments 
following Eq. (51). 

b) The maximum number Nmax of particles below the cut-off is decreasing with in- 
creasing /i', as expected from the fact that the level (n, /) = (30, 30) is getting lower with 
increasing /i' (see Eq. (24) ). 

c) The rms deviations are very small, given the fact that the relevant average energies 
range up to 10^. 

On the contents of Fig. 2 the following comments apply: 

a) The gradual development of supershells with increasing jj is clearly seen in Figs. 2(a)- 
2(d). While in Fig. 2(a) (/x' = 0.0001) no supershcll structure is seen up to = 10000, in 
Fig. 2(b) (//' = 0.001) the first supershell is seen to be completed around = 10000, in 
Fig. 2(c) (//' = 0.002) two supershells are completed up to = 10000, and in Fig. 2(d) 
{jj! — 0.003) most of the third supershell is also completed hy N — 10000. 

b) In Figs. 2(a)-2(g) (/i' = 0.0001 - 0.01) it is clear that A^,i, the value of at which the 
first supershell is completed, is decreasing as a function of /i'. The same is true for Ns2-, the 
value of A^ at which the second supershell is completed. These trends are not followed by 
Fig. 2(h) (//' = 0.02). An explanation of this effect can be gotten from Fig. 1(a), in which 
it is clear that for relatively small jj! {jj! < 0.01, for example) the order of the levels remains 
the same and only their individual energies change, while at fi' — 0.02 and beyond the order 
of the levels changes drastically, especially at higher energies. This drastic mixing of the 
levels at = 0.02 and beyond can also explain the increasing difficulty in determining the 
closing of supershells for large especially at high energies. Fig. 2(h) (/x' = 0.02) being 
the clearest example for this case. 
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c) According to the study by Balian and Bloch of electrons moving in a spherical cavity 
||10|| , which by analogy can be applied to the valence electrons in a metal cluster [Q], the 
minima of the shell energy, Egheiu which correspond to shell closures, should appear at 
equidistant positions (i.e. they should exhibit a periodicity) as a function of N^/^ within 
each supershell. From the data used for plotting Fig. 2 one can see that this condition is 
approximately fulfilled. 

5. Supershells in the 3-dimensional g-deformed harmonic oscillator 

The procedure described in the previous section has been used in exactly the same way 
for the determination of supershells in the case of the 3-dim g-HO. 

The results for Egheii obtained with the 3-dim g-HO for 10 different values of the pa- 
rameter r (assuming ^a; = 1) are shown in Fig. 3, while the relevant parameter values and 
rms deviations are exhibited in Table 2. In order to facilitate comparisons between the two 
models, the first 8 values of r used here are the same as the values of /x' used in the previous 
section. 

As far as the way of calculation is concerned, the maximum shell used in the first 8 
cases was rtmax = 30, as in the previous section, in order to facilitate comparisons between 
the two models. Lower rimax has been used only in the last two shown in Table 2, 

in order to keep the rms deviation small. 

On the contents of Table 2 the following comments can be made: 

a) The behavior of the parameters as functions of r is smooth, even in the case of ag, 
the coefficient of A^^, which in the case of Nilsson's MO was not showing smooth behavior. 

b) The maximum number N^^ax of particles below the cut-off is decreasing with increas- 
ing T. This can be explained by looking at Fig. 1(c). The level used as the cut-off is the 
level with {n,l) = {rimax^^max)^ which, as explained in Sec. 3, is the lowest level within 
the rimax-^^ shell and, as seen in Fig. 1(c), increases very slowly with increasing r. It is 
then clear that the higher the value of r, the more will be the levels coming from below 
and crossing over the cut-off level, resulting in the decrease of Nmax with increasing r seen 
in Table 2. 

c) As in the case of Nilsson's MO, the rms deviations are very small, given the fact that 
the relevant average energies range up to 10^. 
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On the contents of Fig. 3 the following comments apply: 

a) The gradual development of supershells with increasing r is clearly seen in Figs. 
3(a)-3(d) (r = 0.0001 - 0.003), which look very similar to Figs. 2(a)-2(d) (/i' = 0.0001 - 
0.003) of the previous section. 

b) In Figs. 2(a)-2(g) (r = 0.0001 - 0.01) it is clear that A^'^i^g, the value of at which 
the first supershell is completed, is decreasing as a function of r, almost in the same way as 
Nsi is decreasing as a function of fi' in the case of Nilsson's MO, the only difference being 
that for numerically equal values of r and fi' the corresponding value of Ngi^q is slightly 
higher that the relevant value of A^i^, a fact that can be explained by the general property 
of the spectrum of the 3-dim g-HO to expand more rapidly than the spectrum of Nilsson's 
MO, as seen in Sec. 3. The same is true for Ns2,q, A^s3,g, ^^^4,5, i-e. the values of at which 
the second, third, and fourth supershells are completed. These trends are not followed by 
Figs. 3(h) (r = 0.02), 3(i) (r = 0.038), 3(j) (r = 0.05), which correspond to higher values 
of r, where the mixing of the levels is very strong in comparison to the picture existing at 
low values of r (i.e. for r < 0.01, as seen in Fig. 1(c) ). 

c) Despite the fact that the systematics of supershell closures are modified beyond 
r = 0.02 because of the strong mixing of levels, supershells are still seen beyond this point, 
while this was not possible in the case of Nilsson's MO, because of truncation related 
problems, as we have seen in Sec. 3. It should be noticed that this is exactly the region of 
r values which has been found relevant for the description of metal clusters [Q. 

d) In the case of alkali clusters, Na in particular, the first supershell is known to occur 
around = 1000 |]^, [T^, |T^, |T5[, which is quite in agreement with what is seen in Fig. 



3(i), which corresponds to r = 0.038, the parameter value for which the magic numbers of 
alkali clusters are correctly reproduced up to 1500 particles, which is the expected limit 
of validity of theories based on the filling of electronic shells ||^ . It should be noticed that 
only one free parameter exists in the theory, namely r, which has been fixed in Ref. 0] in 
order to reproduce the magic numbers of alkali clusters. Therefore no free parameter has 
been left over in the calculation of supershells. The fact that the present calculation leads 
to a reasonable prediction of the position of the first supershell closure is a quite stringent 
test of the present theory. 
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e) As we have already mentioned, the theory by Bahan and Bloch |]rU| for electrons 
moving in a spherical cavity, which by analogy can be applied to the valence electrons of 
metal clusters predicts that the minima of the shell energy, Esheiu which correspond 
to shell closures, should appear within each supershell at equidistant positions (i.e. they 
should exhibit a periodicity) when plotted versus A^^/^. From the data used for plotting 
Fig. 3 one can see that this condition is approximately fulfilled. 



f) A further prediction of the theory by Balian and Bloch |]T0[ is that the plot of the 

1/3 

cubic root of the magic numbers A^j, corresponding to shell closures, versus the index i, 
counting the shells, should be a straight line In the case of alkali clusters, in particular, 
in which the triangular and squared closed orbits are supposed to dominate |T3|, the 
slope of the line should be 0.603 . In order to make a preliminary estimate of the degree to 
which this prediction is fulfilled, one can employ the data used in plotting Fig. 3(i), which 
corresponds to r = 0.038, i.e. to the parameter value found appropriate for reproducing 
the magic numbers of several alkali clusters. One can then see that this requirement is 
roughly fulfilled, although the slope appears to be gradually decreasing at large i, an effect 
which might be due to missing corrections mentioned in comment h). 

g) In the case of Al clusters, magic numbers have been determined experimentally in 



detail up to 1200 electrons |]I6[, while additional experimental results up to 2700 electrons 
exist The slope of A^j vs. i in this case is considerably lower (around 0.32) [||], 

indicating that closed orbits other than the triangular and squared ones should be present 
in the Balian and Bloch approach ||T^. It will be interesting to examine if the 3-dim g-HO 
can provide any prediction for supershells in Al clusters, after choosing the value of the r 

1/3 

parameter in order to reproduce the right slope of A^ vs. i. 

h) In order to guarantee the reliability of such predictions, it is of interest to study 



in advance the influence of any modifications imposed by the Strutinsky method |Tl|, |12 



which can be applied in cases in which the last shell is open, which are beyond the realm 
of the present analytic study of Section 3. The influence of the well known quantum 
mechanical effect that Uujq should be decreasing with increasing number of particles in the 



cluster [& 18, 191 should also be taken into account. 
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6. Discussion 

The main results of the present study are in summary the following: 

a) The 3-dimensional g-deformed harmonic oscillator, which is known |^ to describe 
very well the magic numbers of alkali clusters up to 1500 particles (the expected limit of 
validity for theories based on the filling of electronic shells), is found in the present study 
able to produce supershell structures, succesfully predicting the first supershell in alkali 
clusters without involving any free parameter in addition to the deformation parameter r, 
which has been fixed for reproducing the magic numbers. 

b) It should be noticed that these successes of the 3-dim g-HO are largely due to terms 
in the Hamiltonian induced by the symmetry, which make the spectrum of the 3-dim g-HO 
to expand with increasing shell number n more rapidly than the corresponding spectrum 
of Nilsson's Modified Oscillator, allowing, among other things, for reliable truncations to 
be performed. 

c) The successful prediction of the magic numbers can be considered as evidence that 
the 3-dimensional g-deformed harmonic oscillator owns a symmetry [the Ug(3) D sOg(3) 
symmetry, which is a nonlinear deformation of the u(3) symmetry of the spherical (3- 
dimensional isotropic) harmonic oscillator] appropriate for the description of the physical 
systems under study. The use of this symmetry for predicting supershells in other kinds 
of metal clusters (Al clusters, for example) is an interesting open problem. The influence 
of using open shells (taken into account by the Strutinsky method [jlT] , p^), as well as of 
the well known quantum mechanical fact that Kuq should be decreasing with increasing 



number of particles in the cluster g, |T^, |I9| , should be taken into account before any final 
predictions can be made. 
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Table 1 



Parameters used for fitting the average part of the total energy (see Eq. (51)) in the case 
of Nilsson's Modified Oscillator, for various values of the model parameter fi' (see Eq. (3)), 
corresponding to the cases exhibited in Fig. 2. The parameters are dimensionless, since we 
have assumed hw — 1 (see Eq. (3)) throughout. The highest shell rimax and the number 
of particles N^ax included in each calculation, as well as the relevant rms deviation cr, are 
also shown. See Section 4 for further discussion. 





"1 


"2 




«4 


10''V/r, 


10^,0 


nm.v 


^ Max 


a 


0.0001 


13.349 


-6.556 


-0.218 


0.972 


4.451 


6.853 


30 


10027 


22.06 


0.001 


-0.025 


0.281 


-1.482 


1.080 


0.080 


-0.132 


30 


10027 


7.12 


0.002 


-0.022 


0.181 


-1.439 


1.073 


0.489 


1.053 


30 


9609 


4.19 


0.003 


-1.281 


1.080 


-1.651 


1.094 


-0.445 


0.282 


30 


9246 


3.40 


0.005 


-1.073 


1.013 


-1.659 


1.096 


-0.623 


0.043 


30 


8641 


2.76 


0.007 


-0.830 


0.910 


-1.654 


1.097 


-0.724 


-0.626 


30 


7937 


4.16 


0.01 


-5.448 


3.432 


-2.155 


1.143 


-2.645 


0.383 


30 


7189 


4.88 


0.02 


5.692 


-3.891 


-0.376 


0.938 


9.135 


-38.541 


30 


4890 


3.50 
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Table 2 



Parameters used for fitting the average part of the total energy (sec Eq. (51)) in the 
case of the 3-dimensional g-deformed harmonic oscillator for various values of the model 
parameter r (see Eq. (29), with q — e^), corresponding to the cases exhibited in Fig. 3. The 
parameters are dimensionless, since we have assumed hcuo — 1 (see Eq. (29)) throughout. 
The highest shell Umax and the number of particles N^ax included in each calculation, as 
well as the relevant rms deviation a, are also shown. See Section 5 for further discussion. 



r 




02 


03 




lO^as 


lO^ae 




Nmux 


(7 


0.0001 


13.585 


-6.554 


-0.228 


0.972 


4.464 
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0.001 


0.753 


-0.061 


-1.424 
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-0.272 
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7.58 
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1.458 


-0.496 


-1.318 


1.059 


2.075 


-1.136 


30 
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4.68 


0.003 


-3.358 


2.186 


-1.855 


1.106 


0.615 


2.501 


30 


9389 


3.52 


0.005 


0.667 


0.166 


-1.500 


1.074 


2.791 


2.422 


30 


8795 


3.49 


0.007 


-2.631 


1.965 


-1.867 


1.106 


2.204 


7.580 


30 
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3.88 


0.01 


-3.657 


2.568 


-2.013 


1.120 


2.523 


15.399 


30 


7893 


6.75 


0.02 


-17.409 


11.661 


-4.267 


1.379 


-10.052 


88.547 


30 


7189 


9.03 


0.038 


-55.417 


41.650 


-13.067 


2.615 


-96.648 


482.365 


26 


4648 


17.98 


0.05 


-64.123 


55.364 


-19.160 


3.777 


-202.994 


999.866 


22 


3020 


16.63 
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Figure captions 



Fig. 1 (a) Energy spectrum of Nilssoii"s Modified Oscillator (in units of huj, see Eq. (3)) 
as a function of the (dimensionless) model parameter jj,'. Magic numbers are shown at the 
main gaps, (b) Energy spectrum of the 3-dimensional g-deformed harmonic oscillator (in 
units of hcuo, see Eq. (29)) as a function of the (dimensionless) model parameter r (with 
q — e'^, where r is real), (c) Same as (b), but extended to higher energy levels. 

Fig. 2 Shell part {Egheii) of the total energy (in units of tiw, see Eq. (3)) for Nilsson's 
Modified Oscillator vs. the number of particles N. The values of the (dimensionless) 
parameter are the same as these listed in Table 1, together with the details of the 
calculation. See Section 4 for further discussion. 

Fig. 3 Shell part {Egheii) of the total energy (in units of hcoo, see Eq. (29)) for the 
3-dimensional g-deformed harmonic oscillator vs. the number of particles A^. The values of 
the (dimensionless) parameter r are the same as these listed in Table 2, together with the 
details of the calculation. See Section 5 for further discussion. 
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